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Objective and Motivation 

A recognized side effect of prolonged microgravity (pg) exposure is visual 
impairment and intracranial pressure (VHP) syndrome [1], Though there is limited 
medical understanding of this phenomenon, it is hypothesized that cephalic shift of 
the cerebrospinal fluid (CSF) and blood in pg may be a contributor. Computational 
models can be used to provide insight into the origins of VHP [1, 2], In order to 
further investigate this phenomenon, NASA’s Digital Astronaut Project (DAP), in 
collaboration with some of the world’s leading experts in ocular biomechanics, is 
developing an integrated fluid physics-based computational model of the human 
body. The model is divided into the eye, the central nervous system (CNS), and the 
cardiovascular system (CVS). This poster su mm arizes our current progress on the 
development and testing of an integrated computational model of CVS and CNS 
based on the whole-body lumped-parameter (LP) model presented by Lakin et al. 
(2003) [3], These models will provide unique capabilities to give insight in 
physiological changes that ca n not be directly measured (e.g. changes in ICP), and to 
help answer questions related to the following research knowledge gaps: 

VIIP1: We do not know the etiological mechanisms and contributing risk factors 
for ocular structural and functional changes seen in-flight and postflight. 

CV7: How are fluids redistributed in flight? 

VIIP13: We need to identify preventive and treatment countermeasures (CMs) to 
mitigate changes in ocular structure and function and intracranial pressure during 
spaceflight. 


Methods 

We reproduced the unsteady 16-compartment whole-body LP model presented by 
Lakin et al. [3] (Figure 1). The general structure of the model: 

• Includes blood, cerebrospinal fluid (CSF), tissue, interstitial fluid, pulmonary 
circulation and organs 

• Captures direct flow between compartments, as well as transfer of fluid 
between capillaries and tissue by filtration 

• Includes compliant interactions between adjacent compartments 

• Incorporates functions to allow the inclusion of the lymphatic system and the 
sympathetic nervous system (SNS) functions 

• Includes a series of 1 3 differential equations to describe the pressure dynamics 
of the system in accordance with the laws of conservation: 

flow in - flow out = rate of volume change (1) 

• Formulates the governing equations in a matrix form: 

W » [f] + [z] * [P] = [I?] (2) 

• Incorporates lymphatic, intracranial, and SNS through forcing and regulatory 
mechanisms 


In addition, the following assumptions are applied: 

• All fluids are assumed to be incompressible and isothermal 

• Pressure-driven flows are laminar and governed by the hydrodynamic 
equation: 

<?.,= p -Yr= z «<, p i- p i) o) 

• Fluid filtration between the capillary and the interstitial spaces is governed by 
Starling-Landin equation 

Filtration — K ct ((P c — P t ) — a ct {Tt c — 7r t )) (4) 


• Compliance between compartments depends on the change in pressure 
differences between compartments in the form of: 


dVij _ „ d( Pi -Pj) 

xr- Ci J~ir~ 


( 5 ) 


• Revokes the Kellie-Monroe Doctrine in order to take into account the 
influence of extra-cranial physiology on ICP dynamics 

• Hydrostatic pressure variation can be approximated through an ad hoc 
variation of resistances in the upper and lower portions of the body 

Legend: 

Cy - compliance between the two compartments Ry - lumped resistance 
K ct - filtration coefficient V tJ - deformation volume of the interface 

P c - capillary pressure between adjacent compartments i and j 

P t - spatially averaged pressure in compartment i Zy - lumped fluidity (1/ Ry) 

Py - spatially averaged pressure in compartment j o ct - capillary membrane reflection coefficient, 

P t - interstitial fluid pressure n c - blood plasma colloid osmotic pressure 

Qy - flow from compartment i to j n, - interstitial fluid colloid osmotic pressure 



Figure 1: Diagram of the unsteady 16-compartment LP model presented by Lakin et al. (2003) [3] 
to simulate pressures, volumes, and flows within the CVS and CNS. C 1J5 Ry, and Qy represent 
compliance, resistance, and flow between compartments i and j, respectively. Numbers in brackets 
represent the spatially averaged mean pressure in mmHg for each compartment. 


Verification and Validation 


The initial verification test challenged the model in a time-dependent mode by 
specifying steady-state values that were offset by 10% from Lakin et aids [3] 
mean values to perturb the system. Since the system returned to the baseline 
mean values, our model successfully conserved all conservation parameters. 

In addition, we performed two validation cases as outlined by Lakin et al. [3] 
to examine the dynamic behavior of the model for the situations of: 

• Pulsatile heart flow, and 

• Postural change 
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Results and Discussions 


The model successfully simulated a pulsatile cardiac cycle that was comparable to 
[3], As can be seen in Figure 2, there are minor discrepancies in the period and 
amplitude of our prediction of heart outflow and Lakin et aids function, as digitized 
from the paper. However, when compared with literature values [4], our simulated 
systolic and diastolic pressures in the central and intracranial arteries were within 
physiological limits of 120/80 mmHg and 100/65, respectively, as presented in 
Figure 3. We therefore considered our model’s response to be satisfactory. 



Figure 2: Plot of heart outflow against time for both the digitized function and the functional 
parameters obtained, as given in [3], 



Figure 3: Simulation outputs for intracranial and central arterial pressure showing diastolic 
pressures within physiological limits of 120/80 mmHg and 100/65, respectively. 

Head-down tilt s im ulations performed by La kin et al. [3] predicted intracra ni al 
arterial flow of 100% and 94% of the mean when the SNS is activated and 
deactivated, respectively. Our implementation of their model predicted <1% change 
with the SNS activated, and almost no change when the SNS is deactivated (Figure 
4). Model verification studies did not identify flaws in our implementation of this 
model. Consequently, we attempted to re-examine the original formulation as 
presented by Stevens and Lakin (2006) [5] but were unsuccessful in resolving the 
discrepancy. Since the acute regulatory model of Lakin et al. [3] still poses 
challenges to general implementation and potentially lacks capabilities critical for 
simulation of chronic pg exposure, we have abandoned their SNS model and will 
develop an independent strategy for regulation based on measured fluid 
redistribution due to long-duration pg [6], 



Figure 4: Plot of the percentage of mean flow from intracranial arteries to intracranial capillaries as 
a function of time. With SNS reflexes activated, the flow remains at approximately 99.99% of its 
mean value. In the absence of these reflexes, the curve tends to 99.34% of its prior value. 

Forward Work 


In order to leverage this model as a foundation for an integrated systems model to 

simulate the CNS and CVS for VHP research, we will: 

• Formally incorporate hydrostatic pressure variation into the matrix equation 

• Modify system parameters such as tissue and flow properties, flowrates and 
pressures to reflect the most current VHP research and to define a reasonable 
physiological envelope that encompasses the astronaut corps 

• Analyze the model via sensitivity studies to identify most significant parameters 

• Test the model against and train the model with independent studies in postural 
change, head-down tilt and pg to develop and validate regulatory functions that 
appropriately capture the observed fluid redistribution 

• Ascertain trends and characteristic ranges of output variables for use by the CNS 
and eye models, including fluid redistribution and cranial blood flow in pg 



